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Abstract 

In a previous paper, we gave a new theoretical framework in which three 
neutrino mixing in matter are discussed. Rigorous analytical solutions are 
obtained. In present paper an approximate method is developed for studying 
three flavor neutrino oscillations in matter in the new framework. Using condi- 
tion of Amf <C AmJ, decoupled resonances, which is discovered by Kuo et al. 
previously, are obtained without small angle approximation. Calculation and 
approximation is consistent and rigorous. All the formulae appear in simple 
symmetric form and some physical characteristic are more explicit. Around 
the two resonant points, the mass eigenvalue whose eigenfunction does not 
take part in the oscillation process, can be easily calculated to higher preci- 
sion. We find, the oscillation amplitude is dominated by vacuum mixing and 
is smaller than 1 for smaller A resonance, but it is not in the larger A. This is 
characteristic for the three flavor oscillations, but do not exist in two flavor. 
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I. INTRODUCTION 



Recently new experimental results for neutrino appear rapidly |T|]. More and more the- 
oretical articles are published attempting at fixing their mass and mixing angles 0. Solar 
neutrino and some other experiments relate to the mixing and oscillations of the neutrino in 
matter For studying these data, theoretical methods for solving the propagation and 

oscillation equations of neutrino in matter are required. Many important results have been 
obtained for two flavors. An important fact, discovered more than ten years ago, is that 
there is resonance of the mixing angle in matter due to the interaction of charged current [|J . 
Because of the increase of mixing parameters the situation is more complex in three flavors 
than two. An important advance was gained in the research of MSW oscillation in three 
flavors in 1986. It was found that the three flavor oscillations are approximately decoupled 
into two independent oscillations and resonances between two flavors ||. Many important 
results from the researches of two flavors oscillation can be used in the three flavor case ||. 
Many works have been done and papers are published. They have been used in analyz- 
ing the experimental results, analytically and numerically However a new theoretical 
framework and analytical method may be helpful to this task. We have made a series of 
efforts in this respect. In a previous paper, the author gave a new theoretical framework 
to discuss three neutrino mixing in matter [EL Results are given analytically in explicit, 
symmetric, and very simple form. Now we apply them to the approximate solutions and 
discuss three flavor oscillations in matter. 

In section II, the results obtained in the first paper M and will be used in the paper, are 
quoted. 

In section III, using and only using the condition of Ami ^ a strict approximate 

mathematical method is developed to deal with the eigenvalue equation of mass in matter for 
three neutrinos. By the method, a cubic algebra equation is reduced into two of quadratic, 
in two regions around A ~ Amf and A ~ Am\ respectively. Thus, we get four solutions 
at each region, but only three of them satisfies our condition imposed on eigenvalues used 
for reducing the equation. In this way, decoupled oscillations and resonances can be read 
out directly. In addition, the solution which does not take part in the resonance can be 
calculated to higher precision which will be very useful in derivating propagating equation 
in a followed paper. 

In section IV, a consistent rigorous mathematical derivation is presented using our the- 
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oretical framework || to expose the behaviors of the mixing parameters in matter around 
A ~ Ami an d A ~ Am?, when Ami <C Am 2 ,. Approximately decoupled oscillation and 
resonance behavior are obtained between and around the A ~ Ami, an d between 
and around A ~ Am\ respectively. This phenomena is the partner of the os- 



had been discovered as the decoupled resonances of 8™ around Ami an d Q™ around Am 2 , 
by T.K. Kuo and his collaborator previously ||. However, there are some phenomena 
which are not discussed or proved explicitly before. In the lower resonance, the amplitude, 



U™ 2 + U%2 2 ~ yVi+vi < 1- ^ is determined by the vacuum parameters. However, in 



the higher resonance, yU^f + U^ 2 ~ 1 which is independent on the vacuum mixing pa- 
rameters. There is notable difference between two flavor and three flavor oscillation and 
resonance in matter, even if in the case of decoupled resonances in three flavors. When the 
vacuum mixing angles are not small, the new observation is meaningful. 

In the end, in section V, the results are discussed and the conclusion is given. 



We begin from the general equation of the eigenvalve problem and some of the rigorous 
analytical solutions obtained in a previous paper || . All symbols and definitions used in that 
paper will be used in here without explanation.. The equation dominated the eigenvalues of 
three neutrinos in matter is. 



dilation and resonance between cos# m and sin8 m in two flavors In three flavor, this 




II. THEORETICAL FRAMEWORK 



A 3 + a\ 2 + b\ + c = 



(1) 



where 



a = — 



(A + Ami) 



(2) 



b = -(Ami) 2 + Am l (Vu ~ V22) + Ama 



(Vu + V22) A 



(3) 



and 



c = (Ami) 2 Am 2 , ~ Ami ^ m 2 (Vu — V22) — Am-^33 A 



(4) 



The solutions for u e mixing matrix elements in matter are 
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3 n 2 

en / / 



1 \ u — Aj 



u = 1,2,3 



where iV" is the normalizing constant. 



Vi 



^ VA M — A,- 



u = 1,2,3 



and A,-s are 



Ai = -\ (m| - ml) = -Amj 

A 2 = \ (ml — mf) = Ami 

A 3 = m\-\ (ml + m?) = Ami 

The mass eigenvalues in matter expressed by \ u are 

M u = X u + l (ml 
Remembering that we have proved 



2 i 2 

7Bn 



^1 + Vl + »73 = 1 



it is easy to see 



jjm2 + ^2 + ^ 



e3 



(5) 



(6) 



(7) 



(8) 



(9) 



(10) 



This relation is helpful to our task below. Now we make use of these formulae and the 
condition Ami "C Am 2 to derive approximate solutions in which decoupled resonances 
emerged explicitly. 



III. APPROXIMATE SOLUTIONS OF EIGENVALUES 



A. Dimensionless Equation 



Before we begin work, it is convenient to make the equation dimensionless in favor of 
estimating the quantitative order for each term in the Eq.(l)-Eq.(4). Taking A, Ami an d 
Ami m t ne unit of A, we introduce the following variable transformation. 
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We get 



A 3 + a\ 2 + b\ + c = 



(12) 



where 



a = — 



(1 + Ami) 



(13) 



b=- 



(Am x ) 2 + Am l (r} U - n 22 ) + Am 2 (7711 + 7722) 



(14) 



and 



^1 ^1 ^1 ^i ^1 

c = (AmJ Am 2 - Am x Am 2 (r/n - r? 22 ) - Ato 1 77 33 



(15) 



Three flavor neutrino oscillations can be decoupled approximately into two independent 
oscillations between two flavors and resonances, if the neutrino masses have the hierarchical 
characteristic, or more exactly if Amf Am 2 , where Amf = (m 2 — ml) /2 and Am 2 = 
m\ — {m\ + m\) /2. This fact have been observed by T.K.Kuo and his collaborator previously. 
Now, we derive it in this new theoretical framework. 



First, we search for the approximate solutions around the region A ~ Am\. In present 

case, 



We try to search for the solutions which satisfy A ~ Am\. Then A ~ Am 1 ~ 1. Neglecting 
all the terms which is of order 1 in Eq.(12)-Eq.(15), we obtain the following approximate 
equation, 



B. Solutions around A ~ Am? 



Aroj ~ 1 and Am 2 ^> 1 



(16) 



1. Solution to the small eigenvalues A ~ Ami ■ 



A" - (7711 + V22) A - {Ami) 2 - Am i (Vn - V22) 



(17) 



There are two solutions for it 



(?7ii + V22) T v fan + V22) +4 (Am x ) 2 - Am 1 (rj u - 7722) 



(18) 



Both of them satisfy the condition A ~ Am 1 ~ 1. Thus we obtain two solutions 



Ai 2 — 



A {rju + r] 22 ) t y A 2 {rju + 7722) + 4 [(Am?) 2 - AAm 2 {rju - rj 22 )} 



(19) 



2. Solution to the large eigenvalue A ~ Am 2 .. 



We try to search for the third solution as A ~ Am 2 ,. That is the solution satisfying 



A ~ Am 2 ^> 1. We can neglect the terms whose orders are equal or smaller than Am 2 in 
Eq.(12)-Eq.(15). We get an approximate equation as 



A 2 - (l + Am 2 .) A + Am 2 (7711 + 7722) = 



There are two solutions 

1 



A 



± 



1 + Am 2 ±J(l + Am 2 ) 2 - 4Am 2 {rju + r/ 22 ) 



(20) 



(21) 



but, only the solution A + satisfies A ~ Am 2 which is used for obtaining the reduced equation. 
Therefore we get a maximum solution 



A3 



1 



A + Ami + \A + Ami) , 1 

v > \ {A 



AAAm 2 {rju + r\ 22 ) 



2^2 



+ Ami) 



(22) 



Remaining to the order A, we get 



A 3 w Am\ + Arj, 



33 



(23) 
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C. Solutions around A ~ Am| 



When A ~ Am\ 



Am\ ~ 1 an d Amj < 1 (24) 



1. Larger solution A ~ Am\. 



We try to search for the solutions it satisfy A ~ Am 2 . Neglecting all the terms having 

/ 2\ 2 

the order smaller than lAm 2 ) , from Eq.(12)-Eq.(15), we obtain an approximate equation 

A 2 -(1 + A^2)A + A^2^ii + ^2 2 )A = (25) 
There are two solutions for it 

A T = 1(1 =f ^(l + A^') 2 " 4^(^11 + ^22) (26) 

Both of them satisfy the condition A ~ Am\. Thus we obtain 



A 2 , 3 = 2 ( A + Am D T~\/(A + Ami) 2 - 4AAm| fa n + 7733) (27) 



2. Smaller solution A ~ Am 
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We try to search for the third solution in A ~ Am 1 <C 1. In Eq.(12)-Eq.(15), neglecting 
the terms whose order are equal or smaller (AmJ , we obtain an approximate equation 



fan + ^22) A - Awij (7711 - 7722) = (28) 

Its solution is 

-2 7/11 - 7/22 



A = Am, 111 ,zz (29) 
?7ii + 7722 
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That is 



Ai = Amf 



2V11 ~ V22 



V11 + V22 



(30) 



In fact, in our present method, higher precision can be reached for this solution by neglecting 
only the terms in order ( Am 1 ) . Now we have a twice algebra equation 



1 + Am 2 )A - Am 1 (r] U - r] 22 ) + Am 2 (7711 + 7722) 



A- 

-2\ 2 



2 2 1 2 / 2 \ ^ 

(AmJ 2 - Am 1 (r/n - 77 22 )J Am 2 + [Am-^J r] 33 = 



(31) 



There are two solutions for this equation, but only one of them satisfies A ~ Am 1 which is 
used to reduce the equation. It is the smaller one. The solution is a minimum solution for 
A ~ Am 2. When we discuss propagating equation in a followed paper, this higher precision 
solution is very useful, but we are not necessary write it here explicitly. 



IV. APPROXIMATELY DECOUPLED RESONANCE BEHAVIORS OF U*{, U*g, 

In section II, we have introduced a set of solutions, Eq.(5) for the U™ u from a previous 
paper ||, because we are interested only in the mixing associated with u e in matter. They are 
expressed as functions of the vacuum parameters and potential A. Now we use the equation 
to discuss the approximately decoupled resonant behaviors of U™{, U^, and U^. 



A. Resonant form for Ai and A2 around A ~ Amf. 

When A ~ Ami, we have get three approximate solutions. They are expressed in 
Eq.(19) and Eq.(22). There is a possible resonance between eigenstates of Ai and A2 when 
the quantity in the square root in Eq.(19) takes the minimum value. It can be reached when 
the following condition of A = A\ is satisfied: 

Ai = 2 Ami VU " V22 2 (32) 
fail + V22) 

Introducing 
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pi = A (7711 + 7722) - 2Amj 



rjn - V22 
Vu + V22 



and 5i = Ami 



2771772 
?7ii + V22 



we can rewrite the solutions Eq.(19) as 



A (7711 + 7722) =F Jp} + 48f 



A ~ Ami 



(33) 



(34) 



It is clear that A5f > when 7^ 0. There is a possible resonance between eigenstates of 
Ai and A2 when the condition A = A\ is satisfied. Because A3 ~ Am^ ^> ~ Ami, this 
resonance is decoupled approximately from A 3 . 

Taking the traditional representation and let U = U^U-iUi ||, we have 



7711 - ^22 



c l c 3 S 1 C 3 



cos26 l i 



(7711 + ^ 22 ) 2 (c?c§ + slclf 

Thus, the lower resonance condition, Eq.(32) , becomes 

9 cos26 l 1 
A = Ai = 2Aml 



r 2 



(35) 



(36) 



It is the same as the result obtained by Kuo et al. 



B. The behaviors of U™, and in the neighborhood of A ~ Amf. 



We can use the equations of U™ and and take the approximation as follows 



UTx,2 « AT 



1.2 



? /-r 



77I 



Ai, 2 + Amf Ai )2 - Amf 



Correspondingly, we have 



iV 



1.2 



V2 



Ai j2 + Amfy vAi i2 — Amf 
Substituting the Eq.(38) into the Eq.(37), we obtain 



jjm 
u el,2 



Ai, 2 - Amf 1 ai^V 7 ?! + ^2 



A li2 -Amf| ,K 2 + (r/f + r/l) 5? 



Where 



ai,2 = V^i + % 2 ( Ai, 2 - ^ , > mj 



^1 -^2 A _2 

2 2"; — 2' 

77f + 77I 



(37) 



(38) 



(39) 



(40) 
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Using Eq.(34) of Ai j2 , we get 



"1,2 = \Jvi+V% {^Pi T \^P 2 i+^ 



(41) 



Then 



jjm 
u el,2 



Ai j2 - AmJ a 1>2 



V1+V2 



it 



pi 



|Ai, 2 - Am? I \a lt2 \\} 2 ^ ' yjfi + 45j 
The 77™ is not independent. It is determined by U™ + + = 1. from it, we obtain 



(42) 



TT m ~ 



(43) 



When the resonance condition pi = is satisfied, we have 

rrm 
u el 



X h2 -Aml a h2 V1+V2 



|Ai >2 - Am? I |tti, 2 



(44) 



C. Resonant form for A2 and A3 around A ~ Am|. 

We have obtained three solutions when A ~ AmJ. They are expressed in Eq.(27) and 
Eq.(30). There is a possible resonance between eigenstates of A 2 and A 3 when the quantity 
in the square root takes a minimum value. It is easy to show that the minimum can be 
obtained when the following condition of A = Ah is satisfied: 

A h = Ami (V11 + V22 - V33) (45) 

Let 

p h = A - Amj (Vn + V22 - 7733) and 5 h = Am^V^n + V22V3 (46) 
the A2,3 can be rewrite as 

A 2 ,s = \ [A + AmlT {pl + 46*) (47) 

It is clear that A5\ > when (r/n + ^22) ^33 7^ 0. We obtain a possible resonance. Because 
Ai <C Ami ~ ^2,3, this resonance between eigenstates of A 2 and A 3 is decoupled approxi- 
mately from Ai. 
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When we take the traditional representation and let U = U2U3U1, we have 



771 = C1C3, 7] 2 = S1C3 and r] 3 = s 3 



Then, 



Vu + V22 ~ V33 = c 2 ± c l + s i c 3 - s 3 = cos 2 #3 



The higher resonance condition, Eq.(45), becomes 

A h = Amj cos26> 3 + O (Am? 
It is the same as the result obtained by Kuo et al. at large A resonance. 



(48) 



(49) 



(50) 



D. The behaviors of U£{, and in the neighborhood of A ~ Am 2 ;. 



When A ~ Am 2 ., the higher resonance, Ai <C A 2)3 . We can use the equations of and 
and make the approximation as follows 



»7i , 772 , 



^3 



.^2,3 ^2,3 A2,3 — Am\ 



Correspondingly, we have 



N 2 



A 



\ 2 / \ 2 

1^ + + 



V3 



A 2 , 3 - Am\ 



Substituting the Eq.(52) into the Eq.(51), we obtain 

A 2 ,3 (A 2 ,3 - Am 2 2 ) "2,3 



U e2,3 



|A 2 ,3(A 2 , 3 - Am 2 ) I / a 2 + p h 



Where 



"2,3 = A 2 , 3 - Ami (vl + vl) 



Using Eq.(47) of A 2j 3, we get 



(51) 



(52) 



(53) 



(54) 



(55) 



Then 
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u: 



A 2 ,3 (A 



2,3 



Am? 



"2,3 



e2,3 



IT 



(56) 



|A 2 , 3 (A 2l3 - Am|) I |a 2j3 | ^ 2 y ^ p 2 + A6h 
The f/^ 1 is not independent. It is determined by + + E/^j = 1. From it, we obtain 

U2 « (57) 



When the resonance condition p ft = is satisfied, we have 



U. 



e2,3 



A 2 ,3 (A 2 ,3 - A^) «2,3 

|A 2 , 3 (A 2>3 - Am|) I |a 2)3 | 



(58) 



V. DISCUSSION AND CONCLUSION 

In situation of three neutrinos in matter, the mixing is generally different to two neu- 
trinos. They do not occur in any plane spanned by two of the three eigenvectors | u m ' u ) 
(u — 1, 2, 3). We can not use a simple mixing angle to describe their mixing as done in two 
neutrino situation. However in case of Am\ <C Ami, we have two approximately decoupled 
resonances. In each case the oscillation do occurs between two mass eigenvectors, but there 
is important difference from two neutrinos situation. 

Around A ~ Ami, matter effect make U™{ and variation when A change. The 
variation has oscillation characteristic and resonance. Oscillation occurs between | v m,x ) and 
| i> rn ' 2 ) but it has an amplitude y/U^ + = yrfi + rfi < 1. This amplitude is determined by 
the vacuum mixing parameters. In plane spanned by | z/ m>1 ) and | u 171 ' 2 ), there exist a mixing 
angle Pi which varies with A. The [7™ = a cos Pi and = asm Pi with a = \Jrjl + r\\ < 1. 
In this meaning, we have oscillation and resonance in a plane spanned by | u m ' 1 ) and | u m ' 2 ). 
It constitutes only a part of | v e ). The | z/ m,s ) contribution is nonzero and has a approximately 
constant mixing coefficient \ u m ' 3 ) ~ r) s | z/ m,s ). Noticeably | z/ m > 3 } enter only due to 
vacuum mixing but not matter effect. It is not oscillation with the potential A change. 
When U e3 = i] 3 is not small, the nature is important. It is different with two neutrino 
oscillation and resonance in matter. 

Around A ~ Am 2 , the case is same. However, because of the small mass, the mixing 
contribution of | z/™' 1 ) is negligible to and oscillation. There exist a angle Ph which 
varies with A in a plane spanned by the eigenvector | u m ' 2 ) and | z/ m,s ). The = cos Ph 
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and = sinfih with a amplitude 1. The oscillation and resonance occur when potential A 
change. In this case, U^l ~ whatever the values of rji is taken. Then the oscillation and 
resonance is more as one of two neutrinos. 

We have used a theoretical framework developed in a previous paper which is convenient 
for dealing with the three flavor neutrino oscillations in matter. In condition of AmJ <C Amj, 
approximately decoupled resonance behaviors is studied and discussed in details. New inter- 
esting physical results are obtained. There are important differences in resonant phenomena 
between two flavors and three flavors even in the case of decoupled resonances, in particular 
when r^3 is large. In traditional mixing angle description, there is mixing order problem. In 
some case, it may mask the physical characteristic and lead to misunderstanding. 
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